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Abstract 

We prove that under certain stability and smoothing properties of the semi-groups generated 
by the partial differential equations that we consider, manifolds left invariant by these flows per- 
sist under perturbation. In particular, we extend well known finite-dimensional results to the 
setting of an infinite-dimensional Hilbert manifold with a semi-group that leaves a submanifold 
invariant. We then study the persistence of global unstable manifolds of hyperbolic fixed-points, 
and as an application consider the two-dimensional Navier-Stokes equation under a fully discrete 
approximation. Finally, we apply our theory to the persistence of inertial manifolds for those 
PDEs which possess them. 
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1 Introduction 

We consider a nonlinear partial differential equation which takes the form of an evolution equation 
having infinite-dimensional Riemannian manifold X as its configuration space. We assume that this 
PDE generates a nonlinear semi-group S{t) which leaves a differentiable compact submanifold 
with boundary M either overflowing or inflowing invariant, and examine the persistence of this 
submanifold under small perturbations of S{t) in a neighborhood of M. Specifically, we prove 
that if a semi-group S'^{t) is within a z/-diameter ball about S{t), and if the submanifold M 
is normally hyperbolic, attracting trajectories at an exponential rate, with a rate of attraction in 
the normal direction greater than in the tangential direction, a continuously differentiable invariant 
submanifold with boundary of the same type Mi, exists for S'^{t) and converges to M as tends 
to zero. Our result generalizes the persistence theory for manifolds invariant to finite-dimensional 
vector fields that has been studied by Moser, Sacker, Hirsh, Pugh, & Shub, and Fenichel (see |44|, 
1 49 1, and |]l9|). Unlike the flows generated by differentiable finite-dimensional vector fields, 
however, the semi-group S{t) does not, in general, possess a bounded inverse and hence, for fixed 
t, the map S{t) is not a diffeomorphism. In particular, its inverse is only defined on the invariant 
manifold. Furthermore, the configuration space is not locally compact, so that one must strongly 
rely on the regularity of the map together with the compactness of M to establish important 
local bounds. 

Recently, Bates, Lu, & Zeng (see have generalized the finite-dimensional persistence theory 
of Hirsh, Pugh, & Shub to semiflows in Banach spaces; in our notation, their configuration space X 
is a Banach vector space. They consider compact manifolds M without boundary as having center- 
type flow, and are able to prove the persistence result for both stable-center and the more difficult 
unstable-center case, the latter being quite complicated due to the general lack of backwards-in- 
time fiow. Herein, our main interest is in applying the general persistence theory to global unstable 
manifolds and inertial manifolds associated with nonlinear partial differential equations, so we will 
restrict our attention to the case of stable flow in the inflnite-dimensional Riemannian manifold. 

We remark that our choice of a Hilbert manifold for the configuration space may be attributed 
to the following reasons. First, we must require the existence of a second order vector field or 
spray on X which generates the geodesic fiow. In general, Banach spaces are not separable so that 
an infinite-dimensional manifold modeled on a Banach space may not be paracompact and hence 
may not have global partitions of unity. Since it is essential for us to be able to construct tubular 
neighborhoods about our embedded submanifold M (we describe these in the next section), we 
must require geodesies on X to exist. In the Hilbert manifold setting, the existence of such geodesic 
flows is equivalent to the existence of a Lagrangian vector fleld associated with the kinetic energy 
Lagrangian. The Riemannian structure generates a natural symplectic form on T*X which is only 
weakly nondegenerate, meaning that the associated linear map taking vector flelds to 1-forms is 
generally not surjective. It is, nevertheless, possible to conflrm that the Hilbert manifolds common 
to most applications do indeed possess geodesies, and herein we shall restrict our attention to such 
examples. A classical example of a configuration manifold which does not have linear structure 
arises in fluids applications. If one is interested in perfect incompressible fluid flow governed by the 
Euler equations, then the approriate configuration space is the closure of the volume preserving 



diffeomorphisms under a Hilbert space norm, a Hilbert manifold (see |17]). Our examination of the 
persistence problem for this PDE shall appear in a future work. Another interesting example is the 
Sine-Gordon equation whose fields take values in S^. One could argue that the analysis of this PDE 
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could be trivially done by working in the universal covering space, however, if one is interested in 
the perturbation of invariant manifolds, and wishes to work in the linear space supplied by the 
cover, one would then have to prove that the perturbed structure remains an invariant manifold 
after projecting from the covering space back to the quotient. 

Another reason for choosing a Hilbert space structure is that we may reduce the normal bundle 
of the invariant manifold M to a Hilbert or unitary bundle, wherein the transition maps are unitary. 
This is essential for the local analysis that we perform. 

Recall that an overflowing (inflowing) invariant manifold with boundary M satisfies S{t)u G M 
for any u £ M for alH < (t > 0) . This means that the infinite-dimensional vector field defined by 
the partial differential operator points outward (inward) on the boundary of the manifold. The exis- 
tence of such an overflowing (inflowing) invariant manifold for the perturbed semi-group S"" (t) is 
established by standard contraction mapping arguments. Specifically, after diffeomorphically iden- 
tifying the infinite-dimensional manifold X with the normal bundle over M in a neighborhood V of 
M, we search for in the space of sections of this normal bundle. The invariant manifold is the 
image of the particular section which is the fixed-point for each time-t map S^{t). We remark that 
Fenichel's proof in finite dimensions in its local form essentially works for the infinite-dimensional 
vector bundle setting with minor modifications, especially for showing that the persisting mani- 
fold is Lipschitz. We provide a different proof, however, that the persisting manifold is actually 
continuously differ entiable. 

As an application of our persistence theorem, we show that the unstable manifold of a hyperbolic 
fixed-point of a PDE satisfles the conditions sufficient for persistence so that the approximating 
semi-group possesses a continuously differentiable unstable manifold of its hyperbolic fixed-point 
which uniformly converges to the unstable manifold of the unperturbed system (see [43| for a proof 
using the deformation method for perturbation of the linearized system). We may then deduce, 
following the work of |2^, that the global unstable manifold of the obtained stationary solution, 
defined by evolving the perturbed local unstable manifold forward in time, is lower semi-continuous. 
In case the global attractor is the closure of the unstable manifolds (of overflowing manifolds) such 
as with gradient systems, we find that the attractor is lower semi-continuous as well. 

As the majority of nonlinear PDEs generate semi-groups which have no explicit representation, 
it is of great interest to consider numerical schemes as the close approximations of the semi- 
group. In fact, our motivation for perturbing the semi-group rather than the PDE may primarily 
be attributed to the inability of numerical schemes to approximate infinite-dimensional vector fields 
in a sense. As an application, we show that the unstable manifold of the hyperbolic stationary 
solution to the 2D Navier-Stokes equation persists under a fully discrete numerical approximation. 

Next we apply our theory to the persistence of inertial manifolds for those PDEs that possess 
them. An inertial manifold for a PDE is a smooth finite-dimensional, exponentially attracting, 
positively invariant manifold containing the global attractor (see |21|, p^]). Subsequent works 
extended the existence results to more general equations and provided alternative methods of proof 
(see for example |48] and the references therein). Most existence proofs require a gap condition to 



hold in the spectrum of the linear term. As shown in |4^, this gap condition also insures that the 
inertial manifold is normally hyperbolic. Thus, our persistence theory applies to PDEs satisfying 
this gap condition; however, since our theory is local, we may only conclude that the perturbed 
system possesses an inflowing manifold that is close to the inertial manifold. In particular, the 
inflowing manifold for the perturbed system may not globally exponentially attract all trajectories, 
as often occurs when discretizing dissipative systems. We include cases where the inflowing manifold 
is an inertial manifold for the perturbed system. 

Our applications complement existing results relating the long-time dynamics of PDEs to the 
long-time dynamics of their numerical approximations as the time and space mesh are reflned. In 
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|]13| , p5|| sufficient conditions are found on a Galerkin scheme approximating the Navier-Stokes 
equations (NSE) to infer the existence of a nearby stable stationary solution for the NSE from the 
apparent stability of the time-dependent Galerkin approximate solutions. Stable periodic orbits 
are studied in ||5^, while hyperbolic periodic orbits can be found in Q. Other aspects of the large- 
time behavior of the Galerkin approximation to the exact solutions of the NSE are considered in 
| |28| ; subsequent results examined the behavior of solutions near equilibria of PDEs under numerical 
perturbation (see [|^], [Q, (3^). For more general global attractors, see [|2^ for sufficient conditions 
on upper semi-continuity. For persistence of inertial manifolds under numerical approximation, see 
m, ll], m, H, H], and II. 

We structure the paper in the following way. In Section 2, we prove that if a nonlinear PDE has 
an overflowing invariant manifold that is normally hyperbolic and whose semi-group satisfies certain 
regularity conditions, then a continuously differentiable overflowing invariant manifold persists 
under perturbation of the semi-group. In Section 3, we state the general form of the nonlinear 
partial differential equation and define the domain of the semi-group which it generates. In Section 
4, we prove that the global unstable manifold of a hyperbolic fixed-point (of the PDE described in 
Section 3) persists as a continuously differentiable overflowing invariant manifold. We apply this 
result to the examination of the two-dimensional Navier-Stokes equation in Section 5, and finally, 
in Section 6, we prove that inertial manifolds persist as well. 



2 Construction of the perturbed manifold 

In this section, we generalize Fenichel's proof of the persistence of overflowing invariant manifolds 
to the setting of an infinite dimensional Riemannian manifold. Much of the proof goes through 
essentially unchanged, particularly the portion which shows that the persisting manifold is Lipschitz 
continuous; however, we provide a complete exposition in order to account for the changes required 
in infinite dimensions, as well as to clarify and generalize some of the details in the technical lemmas. 
As for the final result that the persisting manifold is actually C^, we take a different approach than 
Fenichel, and show that the persisting manifold is a limit of a Cauchy sequence. 

2.1 Infinite-dimensional geometry and sufficient conditions for persistence 

Let (X, g) be an infinite-dimensional C Riemannian manifold modeled on a Hilbert space (H, (•,•)), 
where r > 3. We consider a injective map S : X^X which leaves a submanifold with boundary 
M negatively invariant and has an inverse defined on M. In particular, we assume the existence of 
a compact, connected, and embedded C submanifold M = MU dM such that S~^{M) C M. We 
further assume that M C M', where M' is also a C" connected and embedded submanifold of X 
having compact closure, which is the same dimension as M and satisfies S~^{M') C M. This gives 
each X £ M an M'-open neighborhood, and thus allows us to avoid the complications associated 
with charts defined on domains with boundary. We assume that the differential of S~^\m' is only 
defined in the tangent bundle of M' rather than in the entire tangent bundle of X over M' . 

We consider the tangent bundle {T{X),it) restricted to tt~^{M') and use its Riemannian struc- 
ture to induce the splitting 

T{X)\m' =T{M')®N, 

where is the C"""^ normal bundle (any transverse bundle). It is implicit in our compactness 
assumption that M' is finite dimensional, and for concreteness, we fix the dimension to be m. Due 
to the local trivialization of the vector bundles, there exists an open covering consisting of M'-open 
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sets U = {Ui,i £ 2} of M, such that for each Ui S U, we have the vector bundle morphisms 

r5 : T{M')\u,^Ui x P, r/^ : N\u,^Ui x Q, 

where HI = P © Q. In fact, due to the compactness of M, for any such triviahzing open cover U, 
we may choose finite refined subcovers on which we can define charts: 



(W = {Uli = 1, ...,s],ai),3 = 1, J, UlcUlc-'-d Uf-^ C C//"^ C Uf C U/ 

and ai{U-) = Bj{0), the ball of radius j centered about the origin in P. We simply choose for 
each X e M a chart {Wx,crx) such that Wx C Ui e U and Bj{0) C cTx{Wx) and then define 
Ux = a~^{Bj{0)), and extract the finite cover. Thus, for each i and any p £ U- , (7i{p) is an 
m-vector spanned by the first m elements of a fixed ordered basis of H. 

For the purpose of defining small local neighborhoods in the normal bundle, it will be convenient 
to consider the reduction of (A^, ttItv) to the Hilbert group. Recall, that Hilb{Q) is the subgroup 
of Laut{Q) which preserves the inner product, so that a linear automorphism I is in Hilb{Q} 
if and only if 1*1 = Id on Q. Given our Riemannian metric g, we may construct the Hilbert 
trivializations from our vector bundle trivializations {U^ ,tI^). Namely, for each x G and any qi, 

q2 in Q, there exists a positive definite symmetric operator Tix satisfying gixiQi,Q2) = {TixQi,Q2)- 

1 

Then, the Hilbert trivialization is defined for each i by {U^ ,Ti), where Tj = T^^t^. We have defined 
our Hilbert trivializations over the open covers Uf^^C// of M; thus, any continuous map defined 

on some Uf^iUi can be uniformly bounded when restricted to Uf^^C// ^. Often, we will use such 
compact refinements as domains of sections and operators. 

Henceforth, N will represent the reduction of the normal vector bundle to the Hilbert group. 
We define the bounded subset of the normal Hilbert bundle: 

= {v£Nx: gxiv,v) < e^, x£ M'}. 

defines an open neighborhood of the zero section of the normal Hilbert bundle; the set of locally 
Lipschitz sections of this neighborhood A^^ will contain the negatively invariant manifold of the 
perturbed map S'^ upon identification of an open neighborhood of M' with A^"^. This identification 
is given by the existence of a tubular neighborhood diffeomorphism h taking N^, for e sufficiently 
small, onto V, where V is some X-open neighborhood of M'. 

Lemma 2.1 Let X be a manifold, r > 3, that admits a partition of unity and let M be a closed 
submanifold. Then there exists a tubular neighborhood of Ai in X of class C^~^. 



This result is well known and can be found in most texts on infinite-dimensional manifolds (see 
for example). Of course our Hilbert manifold has partitions of unity since the norm is differentiable; 
furthermore, we have a canonical spray whose geodesic field on T{X) defines the exponential map 
exp. In other words, we require the Lagrangian vector field associated with the kinetic energy of 
the metric to exist. In general, for infinite-dimensional Riemannian manifolds, the metric is only 
weakly nondegenerate, but we will restrict attention to model Hilbert spaces which are isomorphic 
to their dual space. The map h is simply the restriction of the exponential to the normal bundle, 
h = exp\N- Because h is at least a diffeomorphism of class C^, there exists a constant C such that 
\\Dh\\ < C and ||-D/i~"^|| < C onV. In case that C > 1, we will use this constant in certain global 
estimates, so that it may disappear in the local versions of these estimates. We take our closed 



submanifold Ai in Lemma 2.1 to be M' and lei S = h o S o h 
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We note that in many of the apphcations that we have in mind, the manifold X is itself a 
Hilbert space, in which case X is identified with T{X) and the geodesic field is trivial. In that 
instance, the tubular neighborhood is obtained from the smooth diffeomorphism h : N^^X defined 
by (x, f ) 1-^ x+v and we may relax our regularity to r > 2. In fact, if we are content with continuous 
vector bundles over M, we may take r > 1. 

Definition 2.2 Denote by P and Q the following projection operators: 

P : T{X)\m'^T{M'), and Q : T{X)\m'^N. 

We will need the following stability condition on our negatively invariant submanifold. 

Condition 2.3 (Stability) Let B{p) = Q o DS{S-'^{p))\n : Ns-i(^p)^Np, and let i?i E (0,1). 

Then C^\\B{p)\\ < for all p e Uf^^U/'K 

This requirement simply states that to first order, the map S decreases lengths of vectors along 
the normal bundle fibers by at least a factor of C'^i}^ . We show that this implies that the map S 
sends elements of the tubular neig hborhood about S'^{M) C M "closer" to M. 

Lemma 2.4 For < e < Sq sufficiently small, S : A^^| j_2 — >A^'^ '^^^'^\u-'~^ some {^2 £ 

S (U- ) i 

(^1,1). 



Proof. By definition of our tubular neighborhood, every p G M has a normal neighborhood Rp 

diffeomorphic to Bs{0) in TpX. Let p E f^/~^, and let v E A^^|5-i(p). By continuity, we may choose 
e small enough such that h{v) E Rg-i(^p-^ and S{h{v)) E Rp. Since exp'^ o 5* o /i is at least a 
map taking ;Se(0) C Tg-i(^p-^X into TpX, it differs from the linear map DSs~i(^p-^v by 0(||t;|p) terms 
using Taylor's theorem, and thus by C(e^). 

By projecting exp~^ o S o h onto N and using Condition 2.3, we may choose e small enough 



so that C~^t?ie + ©(e ) < C~^-&2£, since C > 1. Then, there exists some p' G RpH M such that 
||(ea;p|Ar)~^|p'(S'(/i(t')))|| = \\Q expp^{S{h{v)))\\. Using the compactness of Uf'"^, we may pick an 

eo small enough which holds for any p E U^^"^ . ■ 

We will use the local Hilbert bundle trivializations that we have constructed to associate to 
each section restricted to some a graph over Bj{Q) in P x Q. This will allow us to compare 
elements from distinct fibers without using parallel translation. For each i E {1, s} and p E C// 
and q E N^, define V'i by letting 

Ti{p,q) = {T^\N{p,q),i^i{p,q)) = {p,Ap{q)), (2.1) 

where ipipiq) ^ Q- We denote by 5^ the sections of A^^ restricted to the cover U^~^ , i.e. 5^ = 
r(A^^| s TjJ-s)- Then for all u E 5^, we may correspond to u\rjj-a the map Ui : Sj_3(0) C P^Q 

z — 1 z i 

defined by Uj = ipi o u o . For any two elements v,w E 5*", v and w are e-C" close if 

||f — w\\o := max sup ||wj(a^) — 'iWj(2^)|| < £• 

l<i<^:reBj,3(0) 

We define 

\\Ui{z) - Ui{z')\\ 



Lip{u) = max sup 

i<i<s^,yei3j_3(0) 
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if it exists, and let 5| = {u S 5^ : Lip{u) < 5}. Then, for any u G 5|, ho Im{u) is a Lipschitz 



continuous manifold of dimension m. It is important to note that in any chart overlap C// ^CiU^ ^ 



i,k e 

\\ui{z) - Ui{z')\\ ^ \\uk{z) - Uk{z')\\ 

sup II -| < oo =^ sup II -| < oo. 

2,2'ea,{c//-3)n<7fe(c//-3) If -^11 2,2'e<7,((7/-3)n<7fc(c/^^-3) - ^ II 



Indeed, because of our reduction to the Hilbert group, there exists a C^' Hilbert automorphism 



Oki on any overlap region Uf-'-^ n U^-^ so that 



i{jkilm{u\jjj-a^jjj-a)) = Oki '^i{I'm{u\jjj-3^jjj-z)), 

and since we were careful to shrink our patches, these Hilbert maps are also Lipschitz in the operator 
norm topology. This is a slight generalization of the finite-dimensional setting wherein on requires 
each element of an orthonormal basis to be a C""~^ function over each patch, and then relies on the 
equivalence of norms to conclude that the frame is C^^. 

Finally, let v,w £ be continuously differentiable, and to each such section v associate the 
s-vector of linear operators {Dvi, Dvg) where 

DvieC\Bj-3{0),C{^,Q)) for i = l,...s. 

Then, v and w are e-C^ close if they are e-C^ close and 

\\Dv — DwW = max sup \\Dvi{x) — Dwi{x)\\ < e. 



Definition 2.5 Let S" : X^X be v-C^ close to S on V, and define S" = h'^ oS" oh. 



Corollary 2.6 For < e < £q and ly > both sufficiently small, 



Proof. This immediately follows from Definition 2.5 and Lemma 2.4. ■ 

Following ||l^, we will use the implicit function theorem to show that the images of Lipschitz 
sections of under S'^ remain sections of the bundle. 



Lemma 2.7 For 5, v, and e taken sufficiently small, S'^ : S 



'5~ 



Proof. The fact that S'^iUl^^U/-') C M and S-_\M) C M' C U^^iC//"^ together with Corol- 
lary 2^ ensures that the maps xiP] ^) = '^\n o S'^ o u o S^^{p) are well defined for all u & 
and all p G UU-~'^. 

We must show that these maps are injective. In any patch Uf~^, we define the local represen- 
tative of X to be Xii^j^i^) = o ^) ^) <5)|(/J-i ° i ^^'^ i^ote that Xi(0, 0,(5) is the inclusion 
map of Bj-2{Q) into P. By the implicit function theorem for Lipschitz continuous maps (see |^6|), 
there exists a neighborhood X of Xi(0,0, (5) in the Lipschitz topology such that each f £ X maps 
I3j-2{0) into P injectively and satisfies 



^j_3(0) C /(^j-2(0)) C /(^j-2(0)) C ^j-i(O). 
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To each f £ X we identify the injective map / on Uf ^ and have 



C /([//"') C f{Uf~^) C Uf~\ (2.2) 

By taking e, 6, and ly small enough, Xii^i'^y^) can be made arbitrarily close to the inclusion map 
and thus is an element of X. Since this is true for each i, we see that each x{']£,^,^) is injective 
and that L)Uf~^ C x(UC//~^). It is for this reason that we have defined our sections on UC//~^. 
■ 

In order to show that has an overflowing invariant manifold, we will prove that the map 
S'^ takes 5| into itself (more precisely, S'^ maps the images of elements of 5|) and is in fact a 
contraction. To do so, we will need to consider the partial derivatives of its local representation. 
For each point {p,v) G -/V^|r/J, let Lj : N^\rTj^¥ x Q be defined by 

i i 

Li{p,v) = {ai{p),i;i{p,v)) = {x,y). 
Then L-i((x,0)) is an element of the zero section of . By compactness, for any i £ 



we may choose the constant C in Condition 2.3 large enough so that we have the uniform bounds 



||D(Tj|| < C and \\Da- ^\\ < C on UUf ^, and by definition of the Hilbert group, ||D?/'j(p, u)|| = 
\\D%b-\MP,v))\\ = 1 for all {p,v) G N%uj. 

Condition 2.8 (Smoothing) Let A{p) = DS-^\m'{p) ■ TpM'^Tg-i^p^M' , and let G (0,1). 

Then C^A{p)\\\\B{p)\\ < i9i for all p G UUf^ . 

In the case that X is locally compact, DS is uniformly bounded and D(S'|^) has uniformly 
bounded inverse on bounded neighborhoods of M' . In the more general case, we may choose a 
sufficiently small neighborhood of M' contained in the tubular neighborhood on which the same 
conclusion holds. 

Lemma 2.9 (Invertability) There exists a sufficiently small X-open neighborhood V of M' on 
which DS is uniformly bounded and D{S\m) has uniformly bounded inverse. 

Proof. Since S is and M' is compact, there exists C > such that ||-D5'|| < C on M'. For any 
S > and eachp G M', we may choose e(p) such that ||-DS'|| < C+5 on i3e(p)(p). Set V = Upe-pBsip), 
where P is a finite subset of points on M' provided by compactness, and e = minpgnp£(p)- 

Our conditions ensure that Sljjr is a diffeomorphism, and since D{S\^) is uniformly bounded 

on M', [D{S\jjt)]~^ exists and is uniformly bounded by the inverse function theorem. Thus, by 
the uniform continuity of the spectrum of DS and the compactness of M' , we may shrink V if 
necessary so that [D{S\jjt)]~^ is uniformly bounded on V. ■ 

For each (x, y) G P x Q satisfying L~^{x, y) G N'^\ ^j_^2 ^g_i^ jjj-2 y can define the local graph 

J 2 J 2 

transform of a section of A^^, represented locally as a graph over U^ , into a graph over U/^ by 
(x, y) ^ C^j^ix, y) where G^^^ : P x Q^P x Q and is defined by 

G^fc = (Gr,S G^,') ^ (a,, o vrU o o Lr^, o S^ o Lt'). (2.3) 

Using this notation, we may translate the result of Corollary into the local form 

\\G'^,„Hx,y)\\<e, (2.4) 
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and we may use Definition 2.5 together with Conditions |2.3| and |2.8| to get for v sufficiently small 



\\D2G\t,\x,y)\\ < ^2, \\{DiG'^,\x,y))-'\\\\D2G'^,'{x,y)\\ < ^2, ^2 G (^i,l). (2.5) 

Since M' is negatively invariant, Z)iG?^^(x,0) = 0, so for r/ > 0, we may choose e and small 
enough such that 

\\D,G-,\x,y)\\<l^. (2.6) 

Finally, Lemma |2.9| gives us 

\\DG'^k{x,y)\\ < C, and \\[D^G'^,\x , y)]-' \\ < C y{x,y) G U{N^\j^^^_^^-^^). (2.7) 

for some bounded constant G. The first inequality is valid for a sufficiently small because of 



the closeness of S'^ with S assumed in Definition 2.5. For the second inequality, we use the 



uniform continuity of the spectrum; namely, if the spectrum of DiG^^ {x,y) does not intersect a 



neighborhood of zero as given by Lemma |2.9| , then we may choose u small enough so that the same 
is true for DiG^i}{x,y). 

In the case that a section u G 5| is continuously differentiable, we have that 

||[DGr,^(x,n,(x))]-i||<C. (2.8) 

We make a remark on the constant C. Our canonical spray defines a unique bilinear form which 
in turn gives rise to a unique Riemannian connection. For simplicity of analysis, we have pulled 
back this connection in each local trivialization by isomorphisms which we have bounded by C. 

2.2 Lipschitz negatively invariant manifolds for the perturbed mapping 



The proof of the following lemma relies heavily on the smoothing condition |2.8| in its local form ( |2.5D 
which is essential for overcoming the possibly large bound C on the norm of the local derivatives 
ofCr^in (^). 

Lemma 2.10 For e, 5, and v sufficiently small, : 

Proof. To prove that S'^{S^) C 5|, it is equivalent to use the local representation (|2.3| ), and show 
that for all u G 5| 

\\G'^,\x,u,{x)) - G'^,\x',u,{x'))\\ < 6\\G'^,\x,u,{x)) - G'^.^x' ,u,{x'm (2.9) 
for x,x' £ L{N^ j_g Q_i,rrj-3j- As usual, we will get a lower bound on the right hand side of 



equation ( ^.91) in terms of the difference of x and x' . We have that 

\\Gik^{x,Ui{x)) - G^ik^ix' ,Ui{x'))\\ > ||GJ'fc^(x,Ui(x)) - G5;^(x',Ui(x))|| 

- \\G'^,Hx',u,ix))-G^,Hx',u,ix') 

= ii'7'iii-r2ii- 



With Lemma 2^ in its local form (^.7]), DiG^^} is invertible, and so we may use the implicit 



function theorem together with a Taylor expansion to get a lower bound on 71; namely, since 
x-x' = [D,G^,^]-\G\,\x, Ui{x)) - G\,\x', Ui{x)) + 0{\\x - x'f ), 



9 



we have that 

IITill > \\[D,G'^,']-'\\-'\\x - x'W + Oi\\x - x'f)- (2.10) 

The upper bound on ||72|| is simply 

IIT2II < 5\\D2G'^k^\\\\x-x'\\+0{\\x-x'f). 

We subtract ||72|| from \\Ti\\, factor the term use the bounds in (^.7]), and take x' 

close enough to x so that 0{\\x — x'\\) < 5C\\x — x'\\. Then 

IITill - IIT2II > \\[DiG^,']-^-\l - 26C')\\x - x'\\. (2.11) 

Let < g < 1, and choose 5 < This ensures that the constant multiplying ||x — is strictly 
positive. 

Next, we estimate 

WG'^.'ixM^)) - G'f.'ix' Mx'm < (||Z)iGr,2|| + 6\\D2G'^,'\\)\\x - x'\\ + 0{\\x - x'f). 

Then, we take u small enough so that (|2.6|) holds with t] < ^7^^ and choose ||x — even smaller 
if necessary so that 0{\\x — x'f) < r]\\x — x'\\. We get 

\\G^,\x,u,{x))-G';,\x',u,{x'))\\ < -imDiG'^,']-'\\-'+HD2G^k'\\ WiDiG'^.YT') 

Q 

■\\G'ik{x,Ui{x)) - G'if}{x',ui{x'))\\, 

and using ( ^.51 ) and (|2.7D completes the estimate. Finally, we note, as in [|l^, that since Bj-^{Q) 
is compact and convex, the estimate we have just derived holds in the entire ball. Thus, from our 
definition of 5|, we have shown that S^{u) is Lipschitz. ■ 

Theorem 2.11 Assume Conditions [g.^ and \2.8^ hold, and let e, 5, and v he sufficiently small. 
Then there exist Lipschitz submanifolds My of X which are negatively invariant with respect to the 
perturbed mappings S'^ . Furthermore, My^M as i^— >0 in C". 



Proof. The proof proceeds in the usual way. We first show that 5^ is a contraction on 5| in 
the C^ topology. Then, since 5| is a metric space, closed under C^ convergence, we appeal to 
the contraction mapping theorem to show that has a fixed point G 5|, and define My = 
h{Im{u-\-^)). 

Thus, we will obtain a uniform bound on the distance between the image of two sections along 
any fiber over \JU^~^. Let u, v! S 5|; then, for any p in some Uj!~^, the proof of Lemma 2.7 allows 
us to choose some i and x,x' in crj {Uf~^) such that 

a,{p) = G'^,\x,u,{x)) = G'^,\x',u',{x')). (2.12) 

A simple estimate shows that 

\\G'^,\x, u,{x)) - G^,\x', n',{x'))\\ < \\D^G'^,^ ■ \\x - x'\\ + 0{\\x - x'f) 

+ \\D2G'^k'\\ ■ ||n,(x) - u',{x)\\ + 0{\\udx) - Kix)f) 
+ \\D2G'^k'\\ ■ \\<i^) - <ix')\\ + 0{\\v!,ix) - v!,ix')f) 
< {2ri + ^25)\\x - x'W + ^2\\ui - u'iWo + 0{\\u^ - u'JI^p.lS) 
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Taking x' close enough to x so that Oi^x — < ^Hl^iGJ'^^] ^|| ^||x — x'|| and using ( 2.1C| ), we 
have 

m\>\\\D^G\^\^\\-^\\x-x\ (2.14) 



By IITill = \\G\^^{x' ,v!^{x')) - G'i^\x' ,u,{x))\\ so we get 

h^\D^G\^\'\\-'\\x - x'W < \\D2G'^,'\\\\u',{x') - u[{x)\\ + Oi\\u',{x') - n^(x)f ) 

+ p2GJ'fc^||||n^(x) - u^ix)\\ + 0(||n^(x) - u,{x)f) 
< C{6\\x - x'W + - u\\o) + 0{\\x - x'f) + 0(||n^ - Ui\\l). 

Then, since u,u' £ 5|, (|2.13| ) permits us to choose 5 and v small enough and 0(11^ — u'\\q) < 
(1 — '!92)||wi — u'iWo to get the desired contraction. 

Finally, since locally, and u'^ are the fixed points of G?^ and G^^, respectively, we may use 
the contraction property of these maps (which we have just proven) to show that u'^^vP in G^ and 
hence that M^J^M in C°. ■ 

2.3 negatively invariant manifold for the perturbed mapping 

If ti G 5| is continuously differentiable and p G Ul^-JJ^~^, then Du{p) £ C{TpM' ,T^^p^N) can 
locally be represented in each C//~^ by 

L'ui gC°(Sj_3(0),£(P,Q)) for i = l,...s. 

Hence, we associate to Dm the s-vector {Dui, Dus) G [C°(^j_3(0), £(P, Q))]* which is a complete 
metric space when normed by 

IlL'ull = max sup ||L'Mj(x)||, (2.15) 

having the topology of uniform convergence. We then define the subset 

r^ = {Te [C0(i3j_3(0),£(P,Q))]^ : IITII < 6}. 

We let u'^ be the Lipschitz section in 5| which defines the negatively invariant manifold of the 
perturbed mapping 5^^. The invariance can locally be represented as 

uUG':,\x,uUx))) = G'^,\x,uUx)). (2.16) 

In order to show that u'^ is actually of class C^, we will construct a Cauchy sequence in 5| which 
converges to a fixed point u'^ £ ^l, and then prove that the sequence is in fact Cauchy in C^. In the 
case that u'^ is Frechet differentiable, we may differentiate ( 2.16| ), and find that its s components 
satisfy 

DuUG',,\x,uUx))) = [D,G'^,\x,u^)) + D2G'i,\x,u'^{x))DuUx)] 

■ [D,G'^,\x,u'^)) + D2G':,\x,uUx))Dunx)r' (2.17) 

= [A)^:(D<)]{(Grfc^(x,<(x))}, 

where the superscript u'^ on A^j refers to the section on which the partial derivatives of GJ'j!,^ and 
G^^^ are evaluated. We note that ( p. 17 ) is well defined because of Lemma 2^ in its local form ( p.8[ ). 
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Thus, in the case that u'^ is C^, we have computed the derivative of the local graph transform, 
which maps elements G Lj(A^^|— /„2 J into Bj^i{0). In order to obtain a global 

representation for the derivative of u'^ we must "piece" together our local representations using a 
partition of unity argument. This, in effect, allows us to use our local vector space structure to 
define a global derivative, without having to analyze the spray associated to the true covariant 
derivative. We will need the following: 

Definition 2.12 A partition of unity of class C" on a manifold M consists of an open covering 
{Ui} of M and a family of functions £ C'"(A^,R) satisfying the following conditions: 

1. for all m G M, (t)i{m) > 0, 

2. the support of (pi is contained in Ui, 

3. the covering is locally finite, and 

4- for each point m £ M, ^ 4>i{m) = 1. 

This will be used to unite the s local representations of the derivative into a single operator 
equation. By Urysohn's lemma, for each p in a C"^ manifold Ai and W C M a neighborhood of 
p contained in a coordinate neighborhood of p which is diffeomorphic to an open ball, there exists 
a neighborhood U oi p with U C W and a C" function : >M with < (p{m) < 1 if m S A^, 
(p{m) = 1 if m G [/, and (j){m) = if m G W^. Thus, choose 7^ = 1 on Uf~^ with supp{'yi) C f//~^ 
and since the covering is locally finite, so is {71}; then, simply let (j)i{p) = 7i(p)/ X]fc=i lk{p) 

to satisfy Definition |2.12| with M = UUf~^. From ( |2.1?| ) it is clear that 

s 

DuUx) = Y,MpX{Du'^){x), k = l,...s, (2.18) 

where p £ M satisfies S'^{u'^{p)) = u'^{aj^^{x)). The patched-together derivative displayed in (2.18) 
serves as motivation for our next two lemmas. 

To simplify our notation for a fixed section u, we shall define drg^{^) = DrG^jJ' {^,Ui{S,)) for 
r,s = 1,2. 

Lemma 2.13 Let = {T^, ...,T}) G T\ choose u G 5|, and let he defined by 

s 

Tiix) =Y,Mp)^UTI){x), (2.19) 



i=l 



where for each p £ , x = G^^^ {ai{p),u{ai{p))) . Then for all e, 6, and v sufficiently small, 

r2 G r-^. 

(We note that since ||T^|| < 6, ( |2.190 is well defined because of Lemma 2^ in its local form (^.8|).) 

Proof. We first show that ||r^|| < 6, and to do so, we use the estimates (p.5|)-( |2?7| ). Since the 
norm on in ( 2.15| ) is defined by computing the maximum over i G {1, ...,s}, it is sufficient to 
show that for all i,ke {1, s}, \\Ali{Tl){x)\\ < 6 for ah x G Uf'^. We have 

Ti = KUTl) 

= [dig^ + d2g^Tl] [d,g^ + d2g'Tl] 
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= [di9^ + d2g'Tl] [Id + d2g'THdig')-'] ' (dig') 
so \\Ti\\ < \\d^g^ + d2gXH ' WidigT'W ■ \\id + d2g^n {dig^)-% Also, 

'"hg 



\\dig^ + d2g^Tl\\.\\{dig^)-^\\ < {\\dig^\\ + \\d2g^\\ ■ \\tI\\] h.. ^1mi-i 



< Cii + 'd25, 

and 

\\Id + d2g^ Tl{dig^)-^\\ <l + C^5 + 0{5^) 

since \\d2g^ ■ • < 1 for 5 small enough. Therefore, 

||T||| < + &5) + ^25 + 0(52) ^ < ^ 

for V and hence r] sufficiently small. 

Next, we must show that T? is continuous. The domain of G^,^ is N^\-^j-2 ^ -,,w-2,. Choose 

k e {1, ...,s} such that for ^ e Bj^siO) ando-^^(^) G Uj^~^, ^ = G^,^^{x,Ui{x)) for some z G {1, 
and X G Sj_3(0) and a^\() G 17^ . We may choose such a fc, because of the nesting property of 
the C/^'s and (^). Choose the largest c > so that ;Bc(0 C ;Sj_3(0). By compactness, we may 
choose c independent of i, k or ^. 

Then for any d > 0, we must show that for 

11^ — ^'11 < c taken sufficiently small, (2.20) 

\\Tm-Tm\\ < d, where C = G^,\x' ,u,{x')). 
This result relies on the continuity of T^. Let 

A%x) = dig%x) + dsg%x)Tl{x), s = 1,2. 

Then 

11^1(6 -TKOII < \\A\x)\\ • \\A\x)-' - A^x'r'W + \\A\x)-A\x')\\ • P^(:r')"l 
< p2(x)|| • \\A\xr'\\ ■ \\A\x) - A\x')\\ . \\A\xr'\\ 
+ \\A\xr'\\-\\A\x)-A'ix')\\, 

where A^ ^ exists due to Lemma |2 



Using ( 2.11 ), we see that ( 2.20 ) implies that — can be made arbitrarily small. Since Hvl^H 
and \\A^ \\ are bounded due to (2^), the continuity of drg^ for r, s = 1,2, and T/ give us the 
result. ■ 

Lemma 2.14 (a) The mapping ^ (pi-^ki ^ contraction on for each u G 5| with a unique fixed 
point T", and (h) the map u i— > T" is continuous. 

Proof, (a) Fix n G 5| and let T" and be in . Then 

\\KUTi:)-KUSt)\\ 

= \\[dig^ + d2gX''][dig' + d2g'T,^r' - [dig^ + d2g^S^][dig' + d2g'sn-'\\ 

= \\[dig^ + d2g^T^][dig' + d2g^T^r' {[dig' + d2g'S^] - [dig' + dsff'?;"]} [dig' + d2g'sr]-' 
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+ {[dig^Tn - [dig'' + ds^'^ri} [diQ^ + d2g^S\ 



1 cfl — 1 1 



Using the same type of estimate as in Lemma 2.13, we take v small enough for a sufficiently small 

rj to get 

||A^,(T,")-A^,(Slf)|| < {\\d2g^\\ \\{d,9^)-^\\+5\\d2g'\\ \\{d,g^r^\\) {I + 0{6)}\\Tr - 
< {{}2 + 0{5))\\T^-Sn<l^\\Tr-Sf\\, /x<l, 

when we take 5 small enough for 0{5) < 1 — '02, and fi is the maximum of all the contraction 
coefficients obtained in each chart. 

(b) For r = 1,2, let T"*^ be the fixed point of Yl'Pi^ki- Choose and in 5| and for all 
n G N, let the sequences {T } and {T } be defined by 



=;0 



=0 



where Tu and Ti. are chosen such that T 



:0 =0, 



Oi\\u^-u%) 



Using part (a), we see that 
the sequences are Cauchy and must have limits and T" , respectively. Thus, in order to prove 

that ||T" — T" II = 0(||n — u ||o), we may use induction, and show that if ||r,^. — Tj. \\ = 
Odlii — u ||o), then ||T;j — Tj;.|| remains Odlu — ||o)- The proof of this estimate is similar to 
those already shown above. ■ 



Theorem 2.15 Assume Conditions 2.S and 2.8 hold, and let e, 6, and u he sufficiently small. 
Then the Lipschitz negatively invariant manifolds My of the perturbed mappings S'^ are , and 
TU^'M as i/^O in . 



Proof. Let u'^ be the Lipschitz section supplied by Theorem |2.11| such that ho u'^\jj- = My, and 
let the components of T" satisfy the functional equation = ^ (pi^'ki {T^ )■ We will show that 
is continuously differentiable and that T"" = Du^ . 

We define the sequence {u"} C 5| iteratively by u^''^'^ = S^{u^) so that locally 



ul+\G^^\x,u:{x)) = G^^\x,uUx)). 



u 2 I 



(2.21) 



If "u" G 5| and is of class C^, then u"'^^ G 5| is also using (|2.8| ) together with an implicit 
function theorem argument. Thus, if we choose G 5| such that Du^ G , {u^ = 0, for example), 
then for all n G N, £>n" G [C°(fij_3(0), £(P, Q))]^ By computing the derivative of (|]2l|) 
that 



we see 



i=l 



From the proof of Lemma |2.13| , it is clear that DvJ^ G F'^ for all n G N. Theorem 2.11 shows that 
is a contraction on 5| so {u"} is Cauchy in 5| and has as its unique limit u'^; therefore, to prove 

that {u"} is Cauchy in C^, it suffices to show that {Dn*^} is Cauchy in F*^. 

We argue as in [[^, Lemma 4.1]. From part (a) of Lemma p. 14 , we may choose, for each 

T"" G F'' satisfying 



n, 



Tf = Y,<t^,Ki:{Tf), 



i=l 
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and define Ei<i = supn,m>N\\T^" — T"™||. Since sup„ „j>jv 11''^'^ 
states tliat E^^O as N^oo. Now, for eacli k £ {1, ...,s}, 



-u 



lo^O as N^oG, Lemma 2.14| (b) 



- < - TfW + \\Tf - TfW + ||Tf 

I ru I I II n> n^il iln^ /x il iln^ 



k Ih 



(2.22) 



and 



rpU 



n + 1 . 



< 



rpu" 
^k 



En 



wliere fj, is tlie contraction coefficient supplied by Lemma 2.14| (a). 
Proceeding by induction, we find tliat for n > m > N , 



+ En 



so using ( 2.22| ), for n > m > N , 



T" II + 



En 
1^' 



< 2^'"-^||Z?n^ 



r-"|| + 2-^ + i?^. 
I- fi 



(2.23) 



(2.24) 



Then, for any e > 0, we may clioose large enough so that each of the last two term on the right 
hand side of ( p. 24 ) are less than e/3. We then choose m > N large enough so that the first term 
on the right hand side of ( p. 24 ) is less than e/3. This shows that the sequence {n"} is C^-Cauchy. 
We conclude, using the completeness of C^, that vP'^vy in 5| and Du^^Du^ in as n— >oo, and 
by passing to the limit in ( p. 23 ) we obtain Du^ = T"". 

Finally, the contraction property of part (a) of Lemma 2.14 assures us that WDu" - Dn°||^0 
as z^— >0, and together with Theorem 2.11, we obtain My^M in C^. ■ 

Next, assume our mapping S arises as the time-T map of a nonlinear semi- group 5^. As in 
|T9|, we define the following generalized Lyapunov-type numbers 



v{x) = limsupj^oo 
e(x)=limsup,_^i°^^ 



if u{x) < 1, 



where At{x) = DS-t\M'{x) : T,M' ^Ts_,(^,)M' and Bt{x) = QDStiS.t{x))\N : Ns_,{x)^N,. We 
have that u and 9 are bounded on M. Furthermore, it is proven in |^] that the generalized 
Lyapunov-type numbers are constant on orbits and hence depend only on the backward limit sets 
on M. In addition, they give us uniform estimates on the norms of At{-) and -Bi(-) as was shown 
by Sacker and later Fenichel (see |49|, |19|) in the following uniformity lemma. 

Lemma 2.16 a) // lim^^oo ll^t(^)ll/^* = V a; G M, then 3 a < a and C such that ||-Bt(x)|| < 
Ca* V x G M, t > 0. 

(h) Further, suppose that a < 1 and lim^^oo ||^t(2;)|| ||-St(x)||^ = V x E M. Then 3 b < b and C 
such that \\At{x)\\\\Bt{x)\\^ < C V x G M, t>0. 

(c) If fix) < o < 1 and9{x) < 6 Vx G Jd , thenXmit^oo ||^t(2;)|| = andlim^^^oo ll-4t(x)||||St(x)||^ = 
uniformly on M 

(d) V and 6 attain their suprema on M. 



We can now state a corollary to Theorem 2.15. 
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Corollary 2.17 Assume the generalized Lyapunov numbers of the semi-group St satisfy y{p) < 1 
and 6{p) < 1 for all p E M' , and that St leaves the submanifold M negatively invariant. Further, 
suppose that there is a semi-group St which is v-C^ close to St on V in some bounded time 
interval. Then for u sufficiently small, there is a submanifold M'^ which is left negatively invariant 
to St , and converges to M as z^— >0 in C^. 



Proof. By Lemma 2.1(:, we may choose T sufficiently large so that with S = St, Conditions 2.3 
and 2^ are satisfied. Thus it follows from Theorem p. 11 that the graph transform associated with 
St has a fixed point in 5f which we call u'^. 



Using the same argument as in Lemma 2/7, we see that for small t > 0, h ^ (St {h{Im{u'^)))) H 
is the graph of a section G 5|. Since h{Im{u'^)) is negatively invariant to S^ , and S^ 



uc// 



commutes with 5^, we see that 



We can then repeat the argument for all t. 



3 The PDE 

We consider a PDE that may be expressed as an evolution equation on a separable Hilbert space 
EI in the form 

f^+Au + R{u) = Q, 

u{Q)=uo. ^^-'^ 

We denote the inner product in EI by (•, •) and norm | • p = (•,•). We assume that A is a densely 
defined sectorial linear operator with compact inverse. Thus it is possible to choose C ^ such that 
all eigenvalues of A := A + (^Id have strictly positive real part. For a > we define A'^ = {A~°')~^, 
where 

i-" = — - / f-^e-'^'dt . (3.2) 
r(a) Jo 



We denote by the domain of (see |2^, [||]). For a = 0, we define A^ = Id. Then P(A") 

is a Hilbert space with the inner product (A^'u, A°'v) and norm \u\a = \A°'u\ for all u,v G T>{A"). 

The operator A generates an analytic semi-group L{t). We assume that the nonlinear term 
R G C^{V{A^),V{A^-^)) and satisfies 

\R{u)\^^p<C + M{p) \/ueV{A"'):\u\^<p 
\\DR{u)\\op = snp\,\^^^\DR{u)v\.,.f3<M(p) Vn G P(i^) : < p , ^"""'^ 

with 7 > and /? G [0, 1) and where C > and M(-) : iR+ is a given monotonically- 

increasing continuous function with M(0) = 0. Furthermore, for all u G 'D{A"') about which 
we linearize R, we require that A + DR(u) be a closed sectorial operator on EI (see for the 
definitions and properties of sectorial operators) with domain ^^{A) which equals 'D{A) since A 
remains a closed densely defined operator on EI. In fact, for 7 G [0, 1] we have 

Lemma 3.1 Let A be a sectorial operator on EI with a (A) > with A^^ bounded, and let B be a 
linear operator on EI with T){A) C T^{B) and which is A-bounded in the sense that for some K > 0, 
\A^~^Bu\ < K\A'^u\ for all u G T^{A) and all (3 G [0, 7). Then A + B is closed in EI with domain 
V{A) =V{A + B). 
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Proof. Let rE„— >x and {A + B)xn^y in H. Then x„ + Bxn^A~^y so that linin^oo Bxn = 
y — Ax. Hence, y — Ax G M and y S IHI so that Ax S EI and thus x E T^{A). Furthermore, x„— >x in 
'D{A) since —Ax^ = Bxn — y and A~^ {limn^^ Bxn — y) = x in M. 

It now suffices to show that Bxn^Bx in EI in order to ensure that the graph oi A + B is closed 
with domain ^[A). We simply estimate 

\Bxn-Bx\ < \\A^-^\\op\A^~^B{xn-x)\ 

< K-\A^Xn-x)\ 

< using a Poincare-type inequality 

< K-\Axn-Ax\, 

and this completes the proof. ■ 



We remark that although it is shown in Henry [|27], proof of Theorem 1.3.2] that A + B has 



a bounded resolvent on the resolvent set, it is not shown that A + B is closed. Therefore, in 



conjunction with Lemma |3.1| , we have that ^4 + i? is sectorial. In the case that the spectrum of A 
does not have strictly positive real part, we analyze A defined as above. 

Finally, we assume that the solution operator of (|3.l| ), denoted by S{t) : 'D(A"')^'D{A'^), for 
t > 0. In particular, S{t) is a strongly continuous semi-group, such that for each t > 0, S{t) is 
continuously differentiable on 'D(A"'), and its Frechet derivative is bounded on any bounded set. 

Partial differential equations that satisfy the assumptions of this section include the two- 
dimensional Navier-Stokes equations, the Kuramoto-Sivashinsky equation, the complex Ginzburg- 
Landau equations, certain reaction-diffusion equations, as well as many other PDEs. 

Perturbation of the PDE. 

We consider perturbations of the mapping G : 'D(A'^)^'D(A'^) obtained as the time-r map 
of the semi-group of ( |3.1D , G{u) := S{t)u for some fixed r > 0. We denote the perturbed mapping 
by G/j and require that for any R > 0, there exist a K = K{R) > such that 

Assumptions Gh 

. - Gh{v)\.y = \E{v)\^ < K{R)h G ^^(0), 

• \\DG{v) - DGh{v)\\op = \\DE{v)\\op < K{R)h Vt; G ^^(0), where DE is the Frechet differen- 
tial of E. 

We remark that Gh may arise from a fully discrete numerical approximation of the PDE. For 
example, let X^'"^ be a subspace of T){A'^). In the case of a finite element-type approximation, X^''^ 
is a finite-dimensional subspace spanned by a family of polynomials. Let be the projection taking 
T>{A'^) onto X^'"' , and let Gh be the finite element approximation to the map G. In order to define 
our approximation on the same space, we define Gh = Gh ° Ph- Thus, Gh ■ 'D{A'^)^'D{A^), and 
for the largest element diameter h taken sufficiently small, satisfies Assumptions Gh by standard 
approximation arguments. 

4 Unstable Manifolds of hyperbolic fixed points 
4.1 Preliminaries 



Assume that u is a hyperbolic fixed point of the PDE (3.1). We first prove the existence of 



a 



local unstable manifold of u which by definition, is an overfiowing invariant manifold of ( p.l| 
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Then, we show that the global unstable manifold of the hyperbolic fixed point, formed by evolving 
the local unstable manifold forward in time, persists under perturbations of the semi-group 
generated by ( |3.lD . 

We define the new variable 

v(t) = u{t) — u . 

For some fixed r > we define the map G{v) = S{t){v + u) — u. Thus, we have translated the 
hyperbolic fixed point to zero. Hence, v{t) satisfies 

^ + Cv + F{v) = 
v{0) = Vo = Uq - u, 

where 

C = A + DR{u), F{v) = R{v + u) - R{u) - DR{u)v , 

and where DR{u)v is the Frechet derivative of i? at u in the direction v. 

Since n is a hyperbolic stationary solution, our assumptions ensure that the operator C is a 
densely defined sectorial operator with a spectral gap about the imaginary axis. We let cti(C) and 
(72{C) be the spectral sets of C associated with the spectrum of C in the left- and right-half planes, 
respectively. Furthermore, we require the cardinality of ai to be finite. We associate the spectral 
projector P : M^M with ai{C) and Q : H^M with <T2 and so we have that M = PM © QM. Using 



the result of Lemma 3.1, we may use standard techniques to verify for the Navier-Stokes equation 

(4.1) 



below, that for any a G (0, 1), there exists r > such that 

||e'^*t;|| < a\\v\\ Vt > r, Vu G PM n V{A'^) 



||e-^*t;|| < a||7;|| Vt > r, G n , 

where || • || = | • I7 for some 7 G (0, 1]. We remark that the estimates ( p. 3D of the nonlinear term 
dictate which 7 to choose for the space 'D{A^). 
Let 

G{v) = Lv + N{v), (4.2) 

where 

Lv = e-^^v, Niv) = [ e~^'^^~'^F{S{s)v)ds , 

Jo 

with S{s) the semi-group of ( |3.1D . It is clear that DN{0) = and that L has no eigenvalues on the 
unit circle. Because of the smoothing properties of the operator exp(— Cr), the map G : 'D(A'^) 

v{A-y). 

We denote the ball of radius p centered about the origin in T>{A'^) by 

B^iO) := {v G P(i^) : \\v\\ < p} . 

In general, our requirements on the nonlinear term and on its linearization are sufficient to ensure 
that G G^{'D{A'^),'D(A'^)) (see Lemma for details) and that positive constants Ki, K2 exist 
such that 

\\nN{v)\\<K,M{p)p yvGB.iO) 
\\nDN{v)\\op := sup||^||=i \\nDN{v)7]\\ < K^Mip) Vv G ^p(O) , ^^-""^ 
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where DN is the Frechet differential of N , M{p) — > as p ^ and TZ is either Id, P, Q. 

Just as in the proof of Corollary ^.17| , we may show that the existence of a local unstable 
manifold for the map implies that the manifold is an unstable manifold for the continuous equations. 
This is quite advantageous when studying the behavior of overflowing sets of ( ^.ij) under numerical 
approximation . 



4.2 Local Overflowing Manifolds for G 

To show that there exists a local unstable manifold of u we consider the linear system Lv = e~'^'^v 
for some fixed r > so that ( |4.1[) holds. By assumption, the space PMn V^A^) is invariant under 
L and PM. r\V{A'^) f\Bp{0) is a C°° overflowing manifold for L for any p > 0. 

Lemma 4.1 Equation ( p.lD has a overflowing invariant manifold W^'^ in neighborhood of u. 



Proof. Matching our notation from Section 2.1, we set our infinite-dimensional phase space X 



H. We set our unperturbed mapping S = L, the linear map, and our perturbed mapping S'^ = G. 
By construction, the subspace PH n V{A^) is invariant to L, so we set M = PH n 'D{A"') n Bp{<d) 
for p > 0. Then M is a negatively invariant (to L) flat manifold of class C°°, and we have the 
splitting TjM = M = T^M ® for ah x G M, where we may identify T^M with PM n 'D{A"i), and 
N with QH n Vi^A'^). Our tubular neighborhood map is simply h{p, v) = p + v. 

We must verify Conditions and with the constant C in those conditions set to 1. Since 
L is linear operator independent of p G PHn 'D{A"'), for Condition we must show that 

\\QLv\\<'di\\v\\ V u G QIHnP(i^). 



This is precisely (O) with -di = a. For Condition |2.8| , we need ||PL~^|pjjp|j,A7ll " IIQ-^II ^ ^i-, which 
is trivially satisfied by (^]^) since < a. 

Thus, we have shown that has an overflowing invariant manifold, and appeal to Corollary 



2.17 to yield the result. 



4.3 Persistence of global unstable manifolds 

Herein, we study the behavior of the set 

W"" := U^eNG'^Ciy^'O (4.4) 

under perturbations of the mapping G. We assume that W"^ is relatively compact and contained 
in Bji{0) for some P > 0. 

We define the asymmetric Hausdorff semi-distance by 

dist(^,P) := sup inf ||a - 6|| , (4.5) 

where A, B are subsets of D(A'^). We recall that dist(^, P) = if and only if ^ C P. 

Theorem 4.2 Let W^'P be the local unstable manifold of the hyperbolic fixed point u of the mapping 
G given by Lemma 4-1, CLnd let Gh be a G^ perturbation of G satisfying Assumptions Gh- Then, 



there exists a G^ manifold Wj^'^ , overflowing invariant to Gh. Moreover, ifW^ is relatively compact 
in P(j4'^), the set 

W^:=UmmG^{Wt^') (4.6) 

satisfies dist{W^, W^) ^ as ^ 0. 
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Remark. The dynamics on may be different than on W}^ and W may be a proper subset of 
W'^ in the hmit /i — > 0. We illustrate this with an example given in the appendix. 

Proof. We first show that the mapping Gh has a negatively invariant manifold. To do so, 



we choose h > sufficiently small and take p in Lemma 4.1 smaller if necessary so that Gh is in a 
small enough G^ neighborhood of L, so that the conclusions of Lemma 4.1 hold for Gh- From this, 
we may conclude that Gh has a unique G^ negatively invariant manifold W^''^. 

Next, we show that Wj^'^^W^'^ as h^O in G^. To do so, we now consider Gh to be a pertur- 
bation of G rather than L. We let X = HI (or D^A^)), the unperturbed mapping S = G, and the 



perturbed mapping S'^ = Gh- We let M = W^'f which by Lemma 41 is negatively invariant with 
respect to G. From, Theorem 2.15 , VF"''' is G^ close to PMn'D{A'^), hence the normal and tangent 



bundles of these two manifolds are G^ close, respectively. This, together with the continuity of 
the norm and the G^ closeness of L and G, ensure the satisfaction of Conditions |2.3| and |2.8| with 



'di = a' , where p is taken smaller if necessary in Lemma 4.1 to satisfy a' £ [a, 1). Thus, applying 



Theorem 2.11, Gh has a negatively invariant manifold W which converges to VF"'^ as h^O in G^. 



With Sg set to the Lipschitz sections over W'"'^, in order to show that W = W^''', we need only 
prove that VF^''' G 5|, and then appeal to the uniqueness of the contraction mapping theorem. Let 
N^'P be the normal bundle of W^'P, and let <^>l'^ : PM n V{A^) n Bp{0)^V{A^) be the G^ section 
identified with the overfiowing invariant manifold W^''' by the tubular neighborhood map. Then, 
each fiber of A^"''' is a submanifold of 'D{A'^). If there does not exist an element of 5| which can be 
identified with Wj^'^, then there is a p G W^'P such that is transversal to N'^'P\p, contradicting 
the fact that Wj^''' is G^ close to W^'P- Thus, W^''' is a section over W^'P. The fact that this 
section is Lipschitz follows from the G^ closeness of W^'P and PM. n V{A'^) n Bp{0)- 

Next, we prove the lower semi-continuity of W^'P. Since is relatively compact in T>(A"'), 
for any e' > 0, there exists € W^, i = 1, ---jN, for some N <oo such that C ufL-^B^^ie' /2). 
Hence, for any u £ W^, Uh £ T^iA^), 

e' 

\\u - UhW < - + \\u' - Uh\\, (4.7) 

for some i, 1 < i < N . 

To proceed with the argument we need to find «^ £ close to in V{A^). Choose n large 
enough so that G~"'Ui £ W^'P for aW 1 < i < N. Due to the global triviality of the unperturbed 
manifold W^'P, we may once again express the sections of the normal bundle of this manifold as 
5= = {^f G C^(M,QE[ n ViA"/) n r(iV^)}. (The argument that we just gave using the underlying 
C°° flat manifold shows that is not empty.) Let the image of ^h in 'S^ be the fixed point of the 
mapping Gh so that Wj^'^ = ho ^/^I^^^tc^, and define u*^ = G^(G~"(u*), <I>/i(G~"('u*))), where h is 
the tubular neighborhood map. 

We obtain from Assumptions Gh that for any w, z £ ^_r(0) 

\\G{w) - Gh{z)\\ < Kl\\w - z\\ + K{R)h , (4.8) 

where we may assume without loss of generality that Kl > 1. If if G W^'P, and z £ 'D{A'^) is 
defined by {Pz,Qz) = (wj^hiw)), then Ww — >0 as h^O since W^'P and Wj^'^ can be made 
arbitrarily G^ close. By assumption [[^^(zu)!! < R for all n > 0, w £ W^'P- With this choice of w 
and z, ( [4.8D implies that ||G/i(2;)|| < R, by shrinking h if necessary. Thus we may continue to apply 
Assumptions Gh- By iterating this procedure we find that 

\\G-{w)-Gl{z)\\ < K2\\w-z\\+KhY,Kl 

1=0 
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< K2\\w-z\\ + (^j^-^jKh. (4.9) 

We apply this last inequality with w = G'^u^ and z G V{A'^) defined by (Pz, Qz) = {w, ^h{w)). As 
before z||^0 as /i— >0. Hence, for h sufficiently small, using ( [4.S| ) we obtain that ||n*— n^|| < e' /2. 
Therefore, ( fl.Tf ) yields 

J 

inf ll-u-'u'^ll < - + < e' . (4.10) 

n''eiy^ 2 

We conclude that dist(VF^, W'') < e' and so dist(VF", W^) ^ as /i ^ 0. ■ 

The previous theorem only makes explicit use of the properties of the overflowing invariant 
manifold, but does not use the properties of hyperbolic fixed points per se; consequently, this result 
generalizes to any overflowing invariant manifold. In particular, we may combine it with a result 
of Hale, Lin & Raugel [p6[ to obtain 



Theorem 4.3 Let A be the global attractor for which we assume to be compact in T>{A 
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Further, suppose that the union of overflowing invariant manifolds which satisfy Conditions 
\2.$ , and their global unstable manifolds, formed by evolving the overflowing manifold forward in 
time, is dense in A. Then, every perturbation o/(3.1), in the sense of Assumptions Gh, has an 



attractor A that satisfies inax[dist{A, A ),dist{A ,A)] as h —> 0. 

Proof. That dist{A,A^) ^ as /i — > follows the proof of the previous theorem. That 
dist{A^ , A) ^ as /i ^ is given in [^]. ■ 

An example of such a system is a gradient system (see [^] ) . 

5 Application to the Navier-Stokes equations. 

In this section we verify that the Navier-Stokes equations (NSE) and their fully discrete approxi- 
mations satisfy the assumptions of Sections |3| and We will consider the two-dimensional NSE for 
a viscous incompressible fluid in a bounded open simply-connected domain O with boundary d^l 
of class C^. We consider no slip boundary conditions. The governing equations are 

If + (n • V)^ = -Vp + uAu + f in 17 X (0, oo) 1 

V • u = in X (0,oo) \ (5.1) 

u\dn = 0, J 

where / = /(x), the volume force, and u > {), the kinematic viscosity, are given. We denote by 
u = u{x, t) the velocity vector, and p = p{x, t) the pressure which are the unknowns. 

5.1 Estimates on the nonlinear term and exponential dichotomies 



Following the notation in |10|, we set 

V = {v.vG (C^)^, divt; = 0}, 

and define 

M = the closure of V in {L^in)f, 

V = the closure of V in {H^{n) f, 
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where L'^{Q) = and H\n), 1 = 1,2,... denotes the usual L^-Sobolev spaces. H is a Hilbert 

space with the inner product and norm. 

Let Pe denote the orthogonal projection in L'^{Q) x L'^{Q) onto H. We denote by A the Stokes 
operator 

Au = —Pm^u, 

and the bilinear operator 

B{u,v)=Pb{{u-V)v) 

for all u,v in V^A) = Vn {H'^{Q) x H'^{Q)). The operator j4 is a self-adjoint positive-definite 
operator with compact inverse. 

We recall some bounds on the bilinear term B[u,v). The first is 

< c|n|^/2|^i/2^|i/2 |^i/2^|i/2|^^|i/2 yu,v,eV{A), (5.2) 
where c will denote an adequate positive constant. Alternatively, we may use Agmon's inequality 

\\u\\L°°{n) < K\u\^/^ \Au\^/^ G V{A), 

to obtain 

\B{u,v)\ < K\u\^/^ \Au\^/^ \A^/^v\ V-u G V{A), v £ V. (5.3) 

Lemma 5.1 Assume that / G V and let S{t) be the semi-group operator of the NSE. Then S{t) 
is strongly continuous and for each t > 0, S{t)u is continuously differentiable on 'D{A) with locally 
bounded Frechet derivative. Furthermore, the nonlinear term of the NSE satisfies l \3.3^ with 7 = 1 
and P = 1/2, and the global attractor is compact in V^A). 

Proof. The Cattabriga-Solonnikov-Vorovich-Yudovich theorem (see [P, [^], |]57|| ), yields con- 
stants Ci,C2 depending on Q such that Ci||n||(j:^i(Q))2 < |A-^/^n| < C2||ti||(/fi(f2))2 for all u G 
p^^i/2^^ giving us the equivalence of norms. In addition, the projection Pjj : {Hq{Q))'^^{H^{Q))'^ 



is bounded (see for example [10|). Hence, we find that for all u,v €V 
\A'/^B{u,v)\^ <K [ V 



duj dvj d^Vi 

^ +Ui ^ 



dxk dxi dxkdxj 



2 

dn , 



and so for all n, t> G T>{A) by a density argument. We use inequalities of the type ( |5.2| ) on the terms 
in the sum and ( ^.3D on the terms Ui gf^g^. to obtain 

\A'/^B{U,V)\ < i^|Al/2n|V2|yln|l/2|^l/2„|l/2|^^|l/2 

+ K|n|l/2 |^^|l/2 |^^|l/2 y^^^ g p(^) _ 

Then, since Ai|n| < \Au\, Xy^\A^^'^u\ < \Au\, where Ai > is the first eigenvalue of A, we obtain 

\A^/^B{u,v)\ < K\Au\ \Av\ yu,v G V{A) , (5.4) 
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and so the inequalities of ( ^^ ) are satisfied with R{u) = B{u,u) — f and 7 = = 1/2. Further, 
for any uq such that \Auq\ < R, we have that |ylM(t)| < K{R) for ah t > and some K{R) > 
(see @, @, and ||). 

The proof that S{t)u is continuously differentiable for all t > is similar to Lemma 14.3 in |10|, 
merely needing modification of the topology from H to T>{A), and also shows that ||^'S'(t)|_g^(o) II — 
M{R) < 00 for all R > 0. Because we have assumed that the force / € 'D{A^/'^), it follows that 
there is an absorbing ball in T>(A^^'^), and hence the global attractor is in fact compact in I'(A) 
(0). ■ 

Let us suppose that u is a stationary solution such that the operator 

Cv = Av + B{u, v) + B{v, u) 

has no eigenvalues with zero real part. As in the previous section we construct the map G as 
follows. We set v{t) = u{t) — u and 

F{v) := B{v + u,v + u) — B{u, u) — DB{u, u)v = B{v, v). 

Then 

G{v) = Lv + N{v) (5.5) 
= S{t){v + u) -u, 



where 



Lv = e-^\ N{v) = r e-^'^^-'^F{v{s))ds. 
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Lemma 5.2 The operator C is sectorial and G : 'D{A)^'D{A). Moreover, there exists a t > 
such that (|4.l| ), (fl.S]) are verified. 



Proof. From Lemma 3T, we obtain that C is closed in EI with domain D{C) = D{A). In order 
to show that C is sectorial we need prove that [C — A)A~'^ is a bounded operator on EI for some 



a E [0, 1) according to [l27|], Corollary 1.4.5]. Using ( p.2|) and (|5.3D , we obtain that 



\{C-A)A 



< \B{u,A-'^/^v)\ + \B{A-'^/'^v,u) 

< c\Au\ \v\ + c\v\ \Au\. 



That C has compact resolvent and hence has a discrete spectrum is given to us by analyticity 



arguments as in [[13|, Lemma 2.1]. Thus, we may appeal to [|27], Theorem 1.4.8] to obtain the 
equivalence of the norms JCf • 1 and \Af^ ■ \ for all P € [0, 1], where we define Ci = C + cid for c> 
such that Re{a{Ci)). 

As in Section 4.1, we associate the spectral projector P with the spectral set ai = Re{a{G)) < 0, 
and the projector Q with the its complement ^2 = cr/ai. Here, o"(C) is the spectrum of G and we 
have the splitting EI = PM QM. 

Using, [|^], Theorems 1.5.3 and 1.5.4], we have that 



\Ae-^^v\ < Cse'^^lAvl \/v G V{A), t > 0, 
\Ae^^p\ < Cge-T^JAp] Vp G PMnV{A), t > 0, 
\Ae-^^q\ < C73t-i/2e-^*l^V2g| Vg G QM n V{A), t > 0, 



(5.6) 
(5.7) 
(5.8) 
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where —a < Re{ai{C)), < 7 < Re{a2) and for some Cj, > 0. We may also write (5. J 



as 



\Ae-'^*q\ < C^t-^/'^e-"<*K\Aq\ for all q G Qnr\V{A) and t > 0, by appealing to the positivity of 
A. Thus given a < 1 there is a r > such that maxjCse"'*"^, C3r~^/^e~'^'^i^} < a and ( [4 .11 ) are 
verified. 

We show that 

\Ae-^*v\ < Kr'^/'^\A^/'^v\ G < t < r. (5.9) 

Indeed, for v S ^(^4) and with p = Pv and q = Qv we have that 

\Cie-^^v\ < \Cie-^*p\ + \Cie-^*q\ 

< (using ( |5.7D and ( |5.8| ) with the equivalence of norms) 

< C3e-"'^K\Cip\ + C3t-^/'^e-"'^K\Cl^'^q\ 

< (since PM D 'D(A) is finite-dimensional) 

< (and for t € (0, r] ) 

< Cst-^/^e-^'K\cl^'p\+Cst-^/^e-^'K\cl/\\ 

< (using the commutivity of C with the projectors P and Q) 

< C^t-^I\-^'K\C\'\\. 

{K is a generic positive constant.) Hence, (|5.9D follows, and we find that 

\AN{v)\ < r ^ \A'/'Fiv{s))\ds < i^T^/' max \Av{t)\\ 

Jo [T - S) ' 0<*<^ 

Standard energy estimates yield an estimate of the form |Ar;(t)| < K{t)\Av{0)\ for < t < r (c./. 



the analysis of [|53|, Theorem 10.2] and ||lO|]) and the first inequality in (|4.3| ) is satisfied since the 
projectors P and Q are bounded on H. 
We also have that 

DF{v)w = DB{v + u)w - DB{u)w = B{v, w) + B{w, v). 

Thus in the same way one finds that 

\\ADN{v)\\op= sup \ADN{v)ri\ < K\Av\ , 
\Av\=i 

and the second inequality in ([4. 3D is satisfied. ■ 



5.2 Perturbations of the NSE 

To illustrate the verification of Assumptions Gh we consider the fully discrete approximation of the 
NSE given by 

u^+i = {Id + AtuA)-'[4, + AtP^if - B{utul))] , 

where G P/vH and where P/v is the projection onto the first N eigenfunctions of A. With T 
chosen so that ([4.1|) holds, we define the map 



{Id + AtiyA^u^ + At ^(Id + AtiyAy-''-^PN[f - B{u]_^,u]^{)] 
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= Gh{4). (5.10) 

where nAt = T. 

Lemma 5.3 The map defined by 

Gh{v) = Gh{PN{v + u))-u 

approximates the map G, given by (|5.5| ), in the norm on any bounded set. In particular, 
Assumptions G^ are satisfied. 

Sketch of Proof. It is sufficient to prove that the maps G{u) = S{T)u and Gh{u) = &{Pnu) 
are close in the norm on i3/j(0) for any i? > and for h sufficiently small (i.e. N sufficient 



large and At sufficiently small). As shown in |33], G and Gh are continuously close. The proof is 
straightforward but tedious. 

To show convergence is not difficult; one copies the proof for the continuous convergence, 
only applied to the linearization of the equations (as above). That is, one shows that 

n 

DGh{u^Q)^J'Q = {Id + Ati/A)-"Pjv/Uo - Y.^Id + AtiyAy-''-^PNDB{PNu'j_^)fi'^_^ 

i=i 

approximates DG{u) in precisely the same way one shows that ( ^.10| ) and S{T) are close. One 
finds that for all u £ ^_r(0) (we omit the details) 

\A{G{u) - Gu{u))\ + \\A{DG{u) - DGh{u))\U < K{R, T)At + ^^^^ . 



As shown in Lemma 5.1 the global attractor is in fact compact in T>{A). Since is part of 



the global attractor, it is relatively compact in ^^{A). Thus we may restate Theorem 4.3 for the 
map Gh{u), (|1^), with h = max{At, A^y^^}: 

Theorem 5.4 Any overflowing invariant manifold satisfying the conditions of Section 2 is captured 
by the fully discrete approximation, ( ^.10| ), h sufficiently small. Moreover, the global unstable 
manifold of this overflowing manifold, formed by evolving the overflowing manifold forward in time 
is lower semi- continuous in the Hausdorff semi-distance sense at h = 0. 

6 Inert ial Manifolds 

In this section we study the persistence of inertial manifolds associated with a nonlinear evolution 
equation in a Hilbert space H, taking the form 

^ + Au + R{u) = 

u{0)=uo, ^ ^ 

under perturbation of the semi-group. The assumptions on A and R{u) are the same as for 
(|3.1|), but in addition, we assume that ( |6.1| ) is dissipative and that A is self-adjoint, meaning that 
u{t) G i3ijy(0) for all t > T*{uo), where i?o > is independent of uq. 

We assume that the eigenvalues of ^, < Ai < A2 . . . < Xj -^00, repeated with their multi- 
plicities, satisfy for any positive K3, K4 

Am > K3, - Am > KiX'^^^ (6.2) 
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for some m > tjiq, tuq > 0. We define Pm to be the projection onto the span of the generahzed 
eigenfunctions of the operator A corresponding to the first m eigenvalues, and set Qm to be the 
projector corresponding to the complementary spectral set. PDEs that are known to satisfy the 
above assumptions include the Kuramoto-Sivashinsky equation, the complex Ginzburg-Landau 
equations, certain reaction-diffusion equations, as well as other PDEs. (See the references that are 
listed in the introduction.) The above assumptions are sufficient to guarantee that ( |6.1D has an 
inertial manifold. Indeed, we recall 

Theorem 6.1 Suppose th at (|6.2D is satisfied. Then for any R > Rq, e > 0, 5 < 1 there exists a 
sufficiently large m so that (|6.1| ) has an inertial manifold representable as the graph of a function 
$ : PmT^{A.'')^QmT^{A'^) contained in Br{0). Moreover, ^ has the properties 

sup ||$(p)||<e, sup \\D^{p)\\op<d. 
\\p\\<R \\p\\<R 



Again, S{t) is the semi-group operator of ( |6J| ) and we set G{u) = S{t)u for some r > 0. As 
usual, we decompose G{u) = Lu + N{u), where 

Lu := e-^^u, N{u) := - [ L{t - s)R{S{s)u)ds . (6.3) 

Jo 



Clearly, the inertial manifold of ( |6.lD is the inertial manifold for G as well, so to study its persistence, 
we consider perturbations Gh of the map G, satisfying Assumptions Gh of Section ^. 

In general, the perturbed mapping Gh are not dissipative even though the map G is assumed 
to be so. Thus although Theorem 2.15 states that Gh possesses an inflowing invariant manifold, 
it may not be an inertial manifold. For this reason, we will restrict our attention to the truncated 
map Gh which we now define. 

Let e : IR-^ ^ [0, 1] be a fixed function such that 9{x) = 1 for < x < 2, e{x) = for x > 4, 
and |0'(a;)| < 2 for x > 0. Define Oji{x) = 6{x/R^). We consider the map 

u = Gh{v):=Lv + N{v)-eR,{\\vf)E{v) 

:= Lv + N^{v), (6.4) 

where L is given (|6.3| ). The map Gh agrees with the map Gh inside the ball 'S^^^^(O) and it is 
dissipative. Moreover, if we require Assumptions Gh to hold on ;B2_Rq(0), then 

\\G{v)-Gh{v)\\<K{R)h, \\DG{v)-DGh{v)\\op<K{R)h "iv^ViA^) . (6.5) 

Unlike the negatively invariant manifolds that we examined in previous sections, an inertial 
manifold is an example of an inflowing or positively invariant manifold. For this reason, we must 



redeflne our deflnitions of the linear operators A and B given in Conditions 2.2 and |2.8| . We set 

A{p) = DG\m'{p) ■■ TpM'^TGip)M' and B{p) = Q o DG{p)\n : Np^Nc^py 



Theorem 6.2 Suppose that the map Gh satisfies (|6.5| ) and the the map G has an inertial manifold 
as given in Theorem 6.1. Then for h sufficiently small, the map Gh has an inertial manifold in 
Bii{Qi) which converges to the true manifold as h tends to zero. 
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Proof. We set X = V^A^)^ our unperturbed mapping S = G|g^^(o)) ^iid o^i^ perturbed mapping 
S'^ = Gh- The unperturbed inflowing manifold M is the inertial manifold for G. In order to 



show that Gh has an inflowing invariant manifold we appeal to Theorem 2.11, and again verify 



Conditions 2.3 and 



Lemmas 5.4 and 5.5 of |45] shows that that the inertial manifold M is 

these lemmas yield 



normally hyperbolic. Keeping the notation in Conditions 2.3 and 2. 

\\A{p)\\ < i^e(^™+-^2^^"'')^ , 



B(p)\\ < i^e-(^™+i-2^2(i+5)A7„-/i) 



for aU p E n Br,,{0). Hence, (|6j) with r 



^m+i verifies Conditions and ^.8 



Thus, Theorem |2.1l| applies and the mapping Gh has an inflowing invariant manifold that 
converges to M in C''. This inflowing invariant manifold is also an inertial manifold for the mapping 
Gh. ■ 

In [33 1 and |34], the perturbed manifold is constructed off of the linear space Pm'D{A'^). As a 
consequence, a restriction on the size of the time step r in the map of equation (|6.3D is required. In 



addition, the dimension of the inertial manifold may be required to increase to insure that certain 
estimates hold. We have avoided both of these issues by constructing the perturbed manifold off 
of the (inflowing invariant) inertial manifold. 



rem 



If the map Gh is dissipative with absorbing set i3ij„(0), then we may take R = 2Rq in Theo- 
6.1 and merely check Assumptions Gh for R = 2Rq. Furthermore, the inertial manifold for Gh 



is an inertial manifold for the map Gh inside Bjig{0). A more generic case is given by 



Lemma 6.3 Suppose that the ball BRf^{0) attracts orbits o/(6T) exponentially with uniform rate. 
Then for any Rh > Rq there exists an Rh > Rf^ > Ro such that for h sufficiently small the 
ball Bri^ (0) is positively invariant under Gh and Bii<^ (0) attracts orbits starting in Br^ (0) of Gh 
exponentially with uniform rate. Moreover, for such h the inertial manifold for Gh is an inertial 
manifold for the map Gh restricted to ;B/j^(0). 



Proof. By assumption we may suppose that \\u{t)\\ <Ro- {R-Ro)e~^^ for all u{t) solving Q 
with ||tio|| < R, and for some A > 0. Assume that Assumptions Gh hold for the given Rh > Rq. 
Let h be taken sufficiently small so that K{Rh)h + Rq + {Rh — i?o)e~^^ < Rh where K{Rh) given 
in Assumptions Gh. Then we have that for u £ Bji{0) with R < Rh 

\\Gh{u)\\ < \\Ghiu) - G{u) + G{u)\\ < K{Rh)h + R^ + {R- R^)e-^- . 

There is a minimum value, R = R'r, for which the above inequality is less than R. It is given by 



Rh 



Ro{l-e-^^) + K{Rh)h 



Of course Rf^^ ^ Rq as h ^ 0. We choose h sufficiently small so that R^ < V^Rq ( the map 
Gh and Gh agree inside B^j^^{0)). Then for any u G Br^^O) the inequality above shows that 
Gh{u) E BR^{d) and G^{u) converges exponentially to Br<^{0) as n 



oo. 



Persistence of inertial manifolds under numerical approximation has been studied elsewhere. 



In [22 1 and [23|, the persistence of an inertial manifold under a Galerkin approximation of the 
underlying PDE is studied, while |15] studies the behavior of the inertial manifold under a time 
discretization. 



In 1 24] and [31[, finite differencing approximations are studied, and in [p3| ] more 
general maps are considered, as well as finite element approximations and their time discretizations. 
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7 Appendix 



We provide an example of the remark we made following Theorem 4.2, in which the sequence of 
unstable manifolds Wj^ do not converge to the true dynamics of in the limit as h tends to zero. 
We consider the vector field expressed in polar coordinates 



(r — l)(rsin^ 1) + hr sin^ 



9 = sin(0) , 

where /i > 0. Notice that for h = there is a stationary solution at r = 1, 6 = n which is not 
hyperbolic. For h > this stationary solution disappears. Notice also that the vector field with 
/i > is a perturbation of the vector field with h = 0. 




Figure 7.1: Lower Semi- continuity of approximate unstable manifold to the true unstable manifold . 

The dashed line in Figure |7[l| represents the flow of the perturbed system (h > 0); the solid 
lines represent the unperturbed system (h = 0). The unstable manifold of the hyperbolic point at 
r = 1, = of the perturbed system now approximates the unstable manifold of the non hyperbolic 
stationary solution of the unperturbed system as well as the unstable manifold of the unperturbed 
hyperbolic fixed point. By preparing the above vector field these manifolds remain bounded. 

Notice that the unperturbed unstable manifold of the hyperbolic point r = 1, 9 = is a proper 
subset of the perturbed unstable manifold in the limit as /i — > 0. 
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